Quantum calculation of dipole excitation in fusion reactions 
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The excitation of the giant dipole resonance induced by fusion reaction is studied with N/Z asym- 
metry in the entrance channel. The Time Dependant Hartree Fock solution exhibits a strong dipole 
vibration which can be associated to a giant vibration along the main axis of the deformed compound 
nucleus. This dipole motion appears to be non linearly coupled to the shape oscillation leading to a 
strong modulation of its frequency. These phenomenons can be detected in the gamma-ray emission 
from hot compound nuclei. 

where h(p) is the mean-field Hamiltonian. We have used 
the code built by P. Bonche and coworkers with an effec- 
tive Skyrme mean- field and SLy^ parameters [23]. 

Ordered collective motions are a general property of 
mesoscopic systems. In metallic clusters, electron vibra- 
tions are plasmon excitations. In atomic nuclei, oscilla- 
tions of protons against neutrons generate giant dipole 
resonances [1,2]. The general way to excite such modes 
is to use rapidly varying electromagnetic fields associated 
with photons or generated by fast electrically charged 
particles. The collective vibrations can also be thermally 
excited as it was clearly demonstrated in the studies of 
the 7 - emission from hot nuclei [3-6]. It has been re- 
cently proposed that fusion reactions with N/Z asym- 
metric nuclei may lead to the excitation of a dipole mode 
because of the presence of a net dipole moment in the 
entrance channel [7-9]. The first experimental indica- 
tions on the possible existence of such new phenomenon 
have been reported in [10] for fusion reactions and in [11] 
for deep inelastic collisions. However, the real nature of 
such a vibration is still unclear both from the experimen- 
tal and the theoretical point of view. In particular only 
semiclassical approaches or schematic models have been 
used to infer the properties of the generated dipole mode. 

In this letter, we present the first quantum calcu- 
lation of pre-equilibrium giant collective vibrations us- 
ing the time dependent Hartree Fock (TDHF) approach 
[12-15,17]. TDHF corresponds to an independent prop- 
agation of each single particle wave function in the mean 
field generated by the ensemble of particles. It does not 
incorporate the dissipation due to two-body interaction 
[19-21], but takes into account one body mechanisms 
such as Landau spreading and evaporation damping [22]. 
The quantal nature of the single particle dynamics is ex- 
plicitely preserved, which is crucial at low energy both 
because of shell effects and of the wave dynamics. More- 
over TDHF is a strongly non linear theory. Hence it can 
exhibit new couplings between collective modes. 

In the time dependent Hartree-Fock (TDHF) ap- 
proach, the evolution of the single particle density matrix 
p{t) — X]^r=i Ifn) {fnl is determined by a Liouville equa- 
tion. 
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FIG. 1. Density plots, projected on the reaction plane, for the 
central collision reaction ^^O +'^^ Mg at 1 MeV/u. Lines corre- 
sponds to equidistant values of the density. 

The effect of the isospin asymmmetry in the entrance 
channel has been first studied in the -|-^° Mg fusion 
reactions at energies close to the Coulomb barrier. Strong 
quantum effects are expected in these mirror-nuclei reac- 
tions leading to the N=Z "^^Ca compound system. The 
density plots obtained in the central collisions at a ki- 
netic energy of IMeV per nuclcon in the center of mass 
frame is presented in Fig. 1. The system fuses and shows 
quadrupole oscillation around a slowly damped prolate 
deformation. Since and ^'^Mg have significantly dif- 
ferent N over Z ratios (respectively 1.5 and 0.67), the 
protons and neutrons centers of mass do not coincide at 
the initial stage of the collision in contrast to what is ob- 
tained in the fusion of two nuclei with the same N over 
Z ratio. 

Let us first start with head on head collisions along the 
X axis. Adequate observables to study the collective mo- 
tion induced in this mechanism are the dipole moment 
Qd and its conjugated quantity Pd- Qd is the net dis- 
tance between protons and neutrons: Qd = ^^{Xp — Xn) 
where Xp and X„ are the proton and neutron's centers of 
mass coordinates. Similarly, Pd = jWzi-^P ~ where 
Pp — pp and Pn = ^he moments of protons 

and neutrons. In Fig. 2 is plotted Qd as a function of 
Pd- As time goes on, we observe a spiral in the collec- 
tive phase space which signals the presence of a damped 
collective vibration. Indeed, it originates from oscilla- 
tions in phase quadrature of the two conjugated dipole 
variables. 

In order to associate the observed vibration with the 
giant dipole resonance (GDR) in the composite nucleus 
we must study the collective vibration of a '^'^Ca nucleus. 
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The GDR associated to the dipole osciUation of a ^°Ca 
ground state has been excited by an isovector dipole field 
and followed using the TDHF approach. The resulting 
Qd moment is plotted as a fimction of time in Fig. 3 — c. 
The period of the observed oscillations is around 80 fm/c 
which corresponds to an energy E = 15.5MeV in good 
agreement with skyrmc RPA calculations [21] and close 
to the experimental value E ~ 20MeV [16]. For the 
^°Ca formed by fusion (Fig. 3 — a) it is around 150 
fm/c. This large difference is due to the deformation 
of the fused system. Indeed, as we can see in Fig. 1 
the compound nucleus relaxes its initial prolate elonga- 
tion along the axis of the collision with a typical time 
much larger than the dipole oscillation's period. During 
the considered time window, the averaged value of the 
observed quadrupole deformation parameter defined by 
e = ^ <Q^°> is around e = 0.23. 




FIG. 2. Time evolution of the dipole vibration. Dipole moment 
Qd and its conjugate Pd are plotted in the phase spa<;e (a) and Pd 
is plotted in function of time (b) . 

A lower mean energy is expected for this longitudinal 
collective motion Eqdrx ^■s compared to the one simu- 
lated in a spherical ^Ca. Following a macroscopic model 
for the dipole oscillation we expect the energy of the GDR 
to evolve with the deformation as 



EgdRx = Egdr (1 - e) 



(2) 



The frequency of the GDR along the elongation axis, 
fuUfill this relation with e « 0.26 in good agreement with 
the observed average deformation. 

The dominant role of the deformation in the lowering 
of the GDR energy can be easily reproduced by perform- 
ing a Hartree Fock calculation of a ^"Ca nucleus in an 
external quadrupole field. This external field is kept dur- 
ing all the dynamical evolution so that the deformation 
does not relax. At time a dipole oscillation is induced 



by a dipole boost. The resulting dipole oscillations have 
a lower oscillation period. Using the deformation param- 
eter e = 0.23 we find that this period is close to the 
one observed in the fusion case and agrees with the phe- 
nomenological relation (2). 
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FIG. 3. Time evolution of the component of the dipole moment 
parallel to the deformation axis (a) and the component perpendic- 
ular in the collision plane (b) for the reaction 20Q_j_20 ^;jjg 
impact parameter 6 = 0/m (solid line) and b = Afm (dashed line), 
(c) GDR in the '"'Ca with a spherical shape 

To get a deeper insight in the dipole oscillation ob- 
served in fusion reactions we have analyzed the time evo- 
lution of its period. From each point on the collective 
trajectory this quantity can be inferred from the time 
needed to reach the opposite side of the observed spi- 
ral. The resulting evolution is plotted in Fig. 4 — a. This 
period shows oscillations too. On the other hand, we cal- 
culated the monopole moment Qoo and the quadrupole 
moment Q20 which are presented in figure 4 — b. We can 
see that those observables oscillations which are almost 
in phase with those in Fig. 4 — a. This points toward 
a possible coupling between the dipole mode and an- 
other mode of vibration. The evolutions of the monopole 
and quadrupole moments are very similar. In particular, 
they have the same oscillation's period around 166 fm/c. 
Therefore, we conclude that they originates from the 
same phenomenon, the vibration of the density around a 
prolate shape. This oscillation modifies the properties of 
the dipole mode in a time dependent way. 

In order to investigate if this density vibration is the 
origin of the observed behavior of the dipole period we 
can model the induced mode coupling. We consider a 
harmonic oscillator with a spring constant which varies 
in time at a frequency w. This is a simple way to take into 
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account the fact that osciUations of the density, which en- 
ters in the mean field potential of the TDHF equation, 
modify the dipole restoring force. This is a way to take 
into account the non linear behaviour of mean field dy- 
namic. 
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FIG. 4. (a) Time behaviour of the dipole period (solid line) 
and its modelisation with the Mathieu's equation (dashed line), 
(b) Time evolution of the monopole (dashed line) and quadrupole 
moments (solid line). 

It is easy to understand that a lower density as well as 
in an elongated shapes induce a weaker restoring force. 
Thus, variations of the density profiles in the TDHF 
equation can be modeled by a corresponding variation of 
the constant of rigidity in a spring Hamiltonian. In such 
a model, the equation of motion becomes the Mathieu's 
equation 

2 

^ + [l+P{^+2 — )cos{iot)]x = (3) 

where ujo is the pulsation without coupling and lu the 
pulsation of the density's oscillation while /? corresponds 
to the magnitude of the induced frequency fluctuations. 
We have computed the numeric solution of equation (3) 
with the typical external frequency of the monopole and 
quadrupole oscillations. The bare frequency ivo and the 
coupling strength (3 have been tuned in order to get the 
same dipole oscillation frequency and typical strength as 
in the full TDHF calculations. Indeed, because of the 
presence of the oscillating term, the observed frequency 
is different from the bare one. The numerical solution 
of the Mathieu equation shows oscillations with an os- 
cillating period which well reproduces our self consistent 
calculations with P = 0.15 and ujo changed by a factor 
1.2 from the observed value i^GDiix (see Fig. 4-a). Prom 



this analysis it appears that the observed dipole motion 
corresponds to a giant vibration along the main axis of 
an oscillating prolate shape. The observed modulation of 
the dipole frequency is a source of additionnal spreading 
of the resonance line shape. 

Using coherent state picture the maximal elongation 
dmax = ~ 1.4/m is related to the number of 

excited phonons by the expression 



no = 



ME, 



GDRx 



2h^ 



(4) 



where Egdrx ^ 8MeV is the energy of the GDR along 
the (a;) axis and M = is the reduced mass of the 

neutron-proton system, m being the nucleon mass. It 
gives no ~ 1.9. Obviously this value is an upper limit be- 
cause the mean field dynamics underestimate the conse- 
quences of the damping since two body dissipation is not 
taken into account. In the same way as ref. [7,10,18,24], 
we can compute the 7-decay probability P-y over its 
mean value in the equilibrium with the expression 



ri + F, 



(5) 



evap 



where Tevap is the rate at which the compound sys- 
tem decay (see ref. [18]), F^ is the decaying rate of 
the phonons and ugdr is the mean number of excited 
phonons in the equilibrium. The latter can be estimated 

by UGDR = 3e where T is the compound nucleus 

temperature. For ^"^Ca at an excitation energy around 
96 McV, we take Fe^^p « 0.76MeV, F^ « 7MeV and 
Ugdr. ~ 0.6. It gives ^ « 1.2 which shows that the 
pre-equilibrium effects are only a correction for this sys- 
tem since the value for an entrance channel quenching 
(no = 0) of the GDR using the isospin symmetric reac- 
tion ^"Ne +20 Ne would lead to ^ « 0.9. It means the 

first chance gamma emission will be slightly enhanced by 
about a factor 20 %. This is in good agreement with the 
values reported in the littcraturc [10]. 

In order to compute the properties of the pre- 
cquilibrium GDR in fusion we shall also study non zero 
impact parameters. In particular the interplay of dipole 
vibration and deformation can be affected by the rota- 
tion. In addition to the center of mass coordinate system 
with X along the beam axis and y perpendicular to the re- 
action plane we defined a new coordinate system x' ,y' , z' 
where x' is the deformation axis and y = y' is the rota- 
tion axis. For the head on head reaction studied before, 
those two frames are the same. In this case, for symme- 
try reasons the dipole moments along the z' = z and y 
axis are zero at every time. 

If we consider now a finite b, the symmetry only forbids 
the vibration to occur along y. In this case, the ampli- 
tude of the oscillations along x' slightly decreases with 
the impact parameter (Fig. 3 — a, dashed line) whereas 
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an oscillation along z' appears (Fig. 3 — 6, dashed line). 
First we can see that the period of the dipole oscillation 
perpendicular to the deformation axis is lower than one 
along the deformation axis. It is another consequence 
for the prolate deformation of the nucleus indeed the en- 
ergy of this excitation along z' should obey the relation 
Egdb.^' ^ Egdr{^ + e)- This relation gives a value of 
the oscillation's period of 67 fm/c which is the observed 
one on (Fig 3 — 6). The oscillation along z' shows a lower 
amplitude than this along x' (approximativcly 1 per cent 
in relative intensity). This oscillation along the z' axis 
results from a weak symmetry breaking due to the rota- 
tion of the system. Indeed the changing to the rotating 
frame x' ,y = y' ,z' leads to the Hamiltonian with a non 
diagonal term 



H' = R{t)HR-\t)+aJy 



(6) 



where H' is the Hamiltonian in the rotating frame, a is 
time derivative of the angle between the two frames and 
Jy is the generator of the rotations around y. The last 
term induces a motion along the z' axis from a dipole 
vibration along x' . This is indeed what is happening 
since initially the dipole mode is along the x' axis. 

If we now compute the number of excited phonons we 
see almost no variation of the total phonon number no 
with b. Therefore the conclusions reached for 6 = can 
be extended to the full range of impact parameters b. 

Let us now study the role of the incident energy. The 
lack of two body damping in TDHF doesn't enable us to 
go up in energy with this system. So we investigate a 
lower energy 0.5 MeV per nucleon which is the coulomb 
barrier. In this case, the number of exciting phonons is 
1.0. This value is lower than one obtained at IMeV/u 
because the maximum dipole moment decreases with de- 
creasing energy. Therefore the number of pre-equilibrium 
phonons is lower at lower energy. 

We have also investigated the role of the mass of the 
reactions partners to study if the pre-equilibrium effects 
survive in heavier compound systems. We have com- 
puted the ^^Ar Ti at 1 MeV per nucleon in the cen- 
ter of mass system (slightly above the hision barrier). 
We have observed exactly the same phenomenology with 
a strong excitation of pre-equilibrium GDR. The calcu- 
lated no = 0.6 is smaller than the hq obtained for lighter 
nuclei reaction. This comes from the fact that the reac- 
tion with small nuclei is less damped than the one involv- 
ing more nucleons. However the pre-equilibrium GDR is 
also strong in the heavier fusion reaction. 

In summary, pre-equilibrium effects related to the N/Z 
asymmetry in the entrance channel have been for the 
first time investigated using the TDHF approach. The 
relation between the mean oscillation's period and the 
shape of the nucleus have been established. Because of 
the deformation the average dipole frequency along the 
deformation axis is lower than the usual GDR one ,while 



we observe a weak vibration at higher frequency in the 
perpendicular direction in the collision plane. In this di- 
rection, the oscillation occurs only for a non zero impact 
parameter. The evolution of the period have been mod- 
eled by a Mathieu's spring model so that the link between 
the observed phenomena and the giant dipole resonance 
have been established. The frequency modulation due 
to this non-linear coupling between modes plays the role 
of an additionnal spreading of the GDR frequency. The 
results shows that pre-equilibrium gamma emission are 
expected in the case of fusion reaction with asymmetric 
N/Z nuclei. 
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stage of this work. Useful discussions with Pieter Van 
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